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, 1. Introduction 



Progress in exactly solving physical models is often tied to advances in the theory 
of orthogonal polynomials and special functions (OPSF). We here present one such 
instance. In looking for spin chains that admit perfect state transfer (PST) [4] - an 
important task in quantum information processing, we have found a new remarkable 
family of orthogonal polynomials defined on a bi-lattice that we have called for reasons 
^ ■ to be explained below, para-Krawtchouk polynomials. 

^ A breakthrough [14] in the theory of orthogonal polynomials (OPs) has been 

realized recently with the discovery of classical OPs that are eigenfunctions of continuous 
or discrete Dunkl operators defined using reflections. These OPs are often referred to 
as -1 polynomials since they arise through g — )■ — 1 limits of g-OPs. Remarkably, this 
has allowed for a complete characterization of the 4-parameter Bannai-Ito polynomials 
in a discrete variable. A significant feature of these classes of OPs is that their spectra 
Xs depend asymmetrically on the parity of the labeling index. This is also observed in 
the spectrum of the parabosonic harmonic oscillator [11] or equivalently in the discrete 
series representation of the related s/_i(2) algebra [15] which has the dual -1 Hahn 
polynomials [16] as Clebsch-Gordan coefficients. 

It is pertinent to signal that Van Assche and collaborators have investigated 
generalizations of Charlier and Meixner polynomials on bi-lattices [9], i.e. lattices 
formed from the union of two regular lattices. They have observed in particular that 
the recurrence coefficients satisfy Painleve equations. 
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A spin chain is said to allow for perfect state transfer (PST) if an excitation at 
one end can be transported to the other end with probability 1 within some time. 
We consider here XX Hamiltonians with nearest-neighbor interactions whose one- 
excitation restrictions are thus Jacobi matrices J [4]. Determining those systems with 
PST amounts to an inverse spectral problem. As is known [4], PST implies that 
the differences between neighboring eigenvalues of J be odd integers (up to a scale 
factor) . The corresponding system is then characterized by the prescribed form that the 
weights of the associated orthogonal polynomials must take. Indeed, the Jacobi matrix 
(recurrence coefficients) and hence the Hamiltonian can be algorithmically found from 
the weights. 

We thus arrived at the para-Krawtchouk polynomials by looking for quantum spin 
chains with PST that have one-excitation energy spectrum of bi-lattice form. 

This paper has (at least) two readerships in mind. The physicists and quantum 
information specialists might focus on sections 2-3 where this new spin chain with PST 
is described. The OPSF experts might want to concentrate on the remaining sections 
where the para-Krawtchouk polynomials are characterized and also put in the context 
of the quadratic Hahn algebra. 

2. Orthogonal polynomials and quantum spin chains with perfect state 
transfer 

Consider the XX spin chain with Hamiltonian 

7V-1 N 

1=0 1=0 
where Ji > are the constants coupling the sites / — 1 and / and Bi are the strengths of 
the magnetic field at the sites I {I — 0,1, . . . , N). The symbols af, af, erf stand for the 
Pauli matrices. 

Introduce the basis vectors 

|e„) = (0,0,...,l,...,0), 71 = 0,1,2,. ..,A^, 

where the only 1 (spin up) occupies the n-th position. In that basis, the restriction J 
of H to the one-excitation subspace is given by the following (iV -|- 1) x (A^ -|- 1) Jacobi 
matrix 

( Bq Ji ... \ 
Ji Bi J2 ... 

\ ... Jm Bn I 

Its action on the basis vectors |e„) reads. 

|en-i) (2.2) 
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Note also that the conditions 

Jo = Jn+i = (2.3) 

are assumed. 

Let Xg, s = 0, 1, . . . , A?" be the eigenvalues of the matrix J. They are all real and 
nondegenerate. Moreover, they are labeled in increasing order, i.e. Xq < Xi < X2 < 

...Xn- 

To the Jacobi matrix J one can associate the monic orthogonal polynomials 
defined by the 3-term recurrence relation 

P„+i(x) + S„P„(x) + UnPn-i{x) = xP„(x), n = 0, 1, . . . , TV, (2.4) 

where C/„ = > and P_i = 0, Pq{x) = 1. 

Pn+i{x) is the characteristic polynomial of the spectral points Xs 

Pjv+i(a;) = (x- Xo){x - Xi) . . . (x - x^)- (2.5) 

The polynomials Pn{x) satisfy the orthogonality relation 

N 

{Xs)Ws = K Snm, (2.6) 

s=0 

where 

K^U,U2...Un. 

The discrete weights Ws > are uniquely determined by the recurrence coefficients 
The perfect state transfer (PST) condition means that there exists a time T such 

that 

e^^-^|eo) = e^^|e^), (2.7) 

where is a real number. In other words, after some time T, the initial state |eo) evolves 
into the state |eAr) (up to inessential phase factor e"^). 

It is well known that the PST property is equivalent to the two conditions [4] : 

(i) the eigenvalues Xs satisfy 

Xs+i -Xs^ ^Ms, (2.8) 

where Mg are positive odd numbers. 

(ii) the matrix J is mirror-symmetric RJR — J, where the matrix R (refiection 
matrix) is 

/ ... 1 \ 
R ... 1 



\ 1 ... / 

Property (ii) is equivalent to the following [17] 
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(ii') the weights Wg (up to a normahzation) are given by the expression 

^ ^ ' 

This means that if the spectrum Xg is given, the weights Wg are then provided by formula 
(2.9) and the mirror-symmetric matrix J can be uniquely restored by a simple algorithm 
[17]. 

3. Finite bi-lattice and para-Krawtchouk polynomials 

Consider a finite bi-lattice of eigenvalues 

= s + ^(7 - 1)(1 - (-1)^-), s = 0, 1, . . . , iV, (3.1) 
where N is odd and the parameter 7 satisfy the restrictions 

< 7 < 2. 

For 7 7^ 1 we have two uniform sublattices that are intermeshed. One sublattice 
corresponds to even s: 

X2s = Is. (3.2) 

The other sublattice corresponds to odd s: 

X2S+1 = 2s 7. (3.3) 

When 7 = 1 we obtain the uniform lattice Xg — s. The PST condition for uniform 
lattice leads to the Krawtchouk polynomials [1]. The orthogonal polynomials Pn{,x) 
corresponding to the weights (2.9) with Xg given by (3.1) will be called para-Krawtchouk 
polynomials. This term is justified by the observation that the spectrum (3.1) coincides 
(when = 00) with the spectrum of the para-bosonic oscillator [11]. 

Direct calculation with the help of formula (2.9) yields two different expressions for 
the weights Wg depending on whether s is even or odd: 

_ 2-^^(1 + 7/2). (-J).(-7/2 - J). 

- {1/2) J s!(l-7/2). ' ^-O'l'---'^ (3.4) 

and 

2-^(1-7/2), (- J).(7/2 - J), 

^^-^^ ^ (1/2), .1(1 + 7/2). ' - = 0,1,---,^, (3.5) 

where J = (A^ — l)/2 is a positive integer (recall that N is odd). 
The weights are normalized in a standard way 

N 

J]w, = l. (3.6) 

s=0 

Equation (3.6) can be verified by observing that 

S'--^ (1/2), H 1-7/2 
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The hypergeometric function in (3.7) can be calculated by using the Chu-Vandermonde 
identity 

which is valid for arbitrary b, c and for positive integers J. Hence 

(1/2) 



E'"^-- ' n/^f' -'A (3.9) 



s=0 



Similarly 



E^^^^i- ''^1^'^' -V2 (3.10) 



and we obtain the desired property (3.6). 

Note that the weights Ws resemble the corresponding weights for the Meixner 
polynomials on bi- lattices [9]. But in contrast to these examples of orthogonal 
polynomials on bi-lattices that have been considered thus far [9] , the para-Krawtchouk 
polynomials have explicit recurrence coefficients J„ and S„: 

7V-1 + 7 , n{N+l-n)((N + l-2ny-^') 

"~ 2 ^ ^n-Jn- 4(7V-2n)(7V-2n + 2) ^ 

There is an even more convenient representation for the recurrence coefficients 

Un = An-lCn, -Bn = ^ Cn, (3.12) 

where 

, (iV-n)(iV-l-2n + 7) n{N+l-2n-j)_ 

2(2n-7V) ' 2(2n-7V) " ^"^-^"^^ 

The para-Krawtchouk polynomials satisfy the difference equation 

E{x)Pn{x+2)+F{x)Pn{x-2)-{E{x)+F{x))Pn{x) = 2n{n-N)Pn{x),{3.U) 

where 

Surprisingly, equation (3.14) puts the para-Krawtchouk polynomials in the category 
of classical OPs on uniform grids [7] . Note however, that the translation x ^ x + 2 
corresponds to the shift Xs — >■ Xs+2- This means that with respect to the grid (3.1) 
equation (3.14) is of 4-th order. 

F{x) vanishes at the two points a; = and x — ^ which corresponds to the smallest 
points of the two sub-lattices. Similarly, E(x) vanishes at the two points ,x = A'" — 1 and 
X — N — 1 + J which coincides with the largest points of the two sub-lattices. 

When 7 = 1 we obtain the classical symmetric Krawtchouk polynomials with 
recurrence coefficients [5] 

N n(N + l-n) 
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Equation (3.14) becomes in this case the square of the eigenvalue equation for the 
Krawtchouk polynomials [5]. 

Returning to the PST property, consider the affine transformation of the spectral 
points (3.1) 

Xs = axs + /3 (3-17) 

with arbitrary real parameters a, (3. The corresponding monic polynomials Pn{x) will 
satisfy the recurrence relation 

Pn+l{x) + BnPnix) + C/„P„-l(x) = xPn{x), (3.18) 

where 

Using appropriate values for the parameters a,/?, we can always achieve condition (i) 
for the PST iff 

7=^, (3.19) 
where Mi,M2 are positive co-prime integers and Mi is odd. 

4. Relation with the complementary Bannai-Ito polynomials 

In order to obtain an explicit expression for the para- Krawtchouk polynomials we notice 
a remarkable relation with the complementary Bannai-Ito polynomials [14] Wn{x) that 
depend on 4 parameters ri,r2, pi, P2 and are defined through the 3-term recurrence 
relation 

Wn+i{x) + (-l)>2W^n(2:) + Vr,Wn_i{x) = xWr^ix), (4.1) 

with 

n{n + pi - ri + l/2)(n + pi - r2 + l/2)(n - ri - ra) 



V2n 



V2n+1 



{2n + l + g){2n + g) 

(4.2) 

{n + g + l){n + pi + P2 + l)(n + P2 - n + l/2)(n + p2 - ^2 + 1/2) 



{2n + l + g){2n + g + 2) 
where we denote 

g = pi + P2 - ri - r2. 
They can be expressed in terms of the Racah polynomials as follows [14] 



and 



w f\ m T? f -n,n + g+l,p2 + x.,p2-x \ 
VPi + P2 + l,P2-ri-hl/2,p2-r2 + l/2 / 

w ^ ^ (2)^ f -n,n + g + 2,p2 + l + x,p2 + l-x ,\ 

W2n+l[x) ^ K\^>{x-p2hF3[ . 1 , 4.4) 

VPi +P2 + 2,p2 + 3/2,p2 -^2 + 3/2 J 

where the normalization coefficients Kn' are needed to ensure that polynomials Wn{x) 
are monic. 
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Put 



r2 = P2 = 0, ri = '-, pi = ^ ^ , (4.5) 



the recurrence relation (4.1) then reads 

Comparing (4.6) with the recurrence relation for the para-Krawtchouk polynomials we 
conclude that 

Pn{x) = 2" Wn{x/2 -{N-1 + 7)/4) (4.7) 

i.e. that the para-Krawtchouk polynomials can be expressed in terms of the 
complementary Bannai-Ito polynomials with a shifted argument. 

Consider the polynomial Pn^i{x). By assumption, = 2M + 1 is odd, hence by 
(4.3) we have 

/-M-l,x/2-M/2-7/4,-x/2 + M/2 + 7/4 \ , 

Pn+i{x) = const 3F2 ( -y-2M 7+2M ; 1 ) •( 

V 4 ' 4 / 

On the one hand, the hypergeometric function 3^2(1) in (4.8) can be simplified using 
the Pfaff-Saalschiitz summation formula [5] . This yields the factorization 

fl -x\ fx- 2M\ , , 

Pn^x{x) = const , (4.9) 



\ ^ / M+l \ ^ / M+1 

where {x)n = x{x = 1) ... (a; + n — 1) is the shifted factorial (Pochhammer symbol). 

On the other hand, P/v+i(a;) is the characteristic polynomial (2.5). Comparing (4.9) 
and (2.5) we arrive at the explicit expression (3.1) for the spectral points. 

5. Algebraic interpretation. Hahn algebra 

It is possible to relate the para-Krawtchouk polynomials to the quadratic Hahn algebra 
[3], [18], [6]. 

Introduce two operators X and Y on the space of polynomials f{x): let X be the 
operator of multiplication by the argument x : 

Xf{x)^xf{x), (5.1) 

and Y be the difference operator: 

Yf{x) = E{x)f{x + 2) + F{x)f{x - 2) - {E{x) + F{x))f{x), (5.2) 

where the functions E{x),F{x) are given by (3.15). Note that the operator Y is the 
operator on Ihs of equation (3.14) satisfied by the para-Krawtchouk polynomials. 
Introduce also a third operator Z which is the commutator of X and Y 

Z = [X, Y]=XY - YX. (5.3) 

It is then easily verified that the commutators [F, Z] and [Z, X] are quadratic expressions 
in terms of the operators X, Y: 

[y, Z] = -4{X, F}+CiX+GF+Mi, [Z, X] = -4X^+C2Y+GX+M2{5A) 
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where {X, Y} stands for the anticommutator and 

Ci = 4(1 - N^), C2 = -4, G = 4(7V + 7 - 1), (5.5) 

Ml = 2(Ar + 7 - 1){N^ - 1), M2 ^ 2{1 - N){N + - 1). 

The polynomial algebra with 3 generators satisfying (5.3), (5.4) as defining relations is 
known as the Hahn algebra [3]. 
It has for Casimir operator 

Q = Z^-A{X'^,Y}+{Ci+lQ)X^+C2Y^+G{X,Y}+{2Mi-AG)X+2M2Y{5.6) 

which commutes with the operators X, Y, Z: [Q, X] = [Q, Y] = [Q, Z] = 0. 
In the reahzation (5.1), (5.2) Q reduces to a constant 

Qf{x)^qf{x), (5.7) 

where 

(N -1)(N + ^ -1)(N^ -2N + -f(N + 3) - 7). (5.8) 
The Hahn algebra admits a basis 7r„, n = 0,1, . . . , N, where the operator Y is diagonal 
YtIu = A„7r„ (5.9) 
whereas the operator X is 3-diagonal 

XtTu — an+l^n+l + BnTTn + an^n-l- (5.10) 

The coefficients a„, S„ and A„ can be calculated from the representations of the Hahn 
algebra [3] with the initial conditions Qq — a^+i — 0: 

Xn = 2n{n- N). (5.11) 

The coefficients and -B„ coincide with the recurrence coefficients (3.11) for the para- 
Krawtchouk polynomials. 

This gives a simple algebraic interpretation of the para-Krawtchouk polynomials. 
It is obvious that the basis vectors 7r„ are given by the para-Krawtchouk polynomials 
in the realization (5.1), (5.2): 

7rn = Pnix). (5.12) 

It is worth noting that the para-Krawtchouk polynomials can be obtained from 
a periodic reduction (with period 4) of the Darboux transformations of finite Jacobi 
matrices. Indeed, as was already noticed in [12], such periodic reductions will lead 
to the Hahn algebra and to the appearance of nontrivial sub-lattices for the spectral 
points. Only ordinary Hahn polynomials were analyzed from this point of view in [12]. 
The para-Krawtchouk polynomials provide a new non-trivial example of such periodic 
closures. 



Para-Krawtchouk polynomials 



9 



6. The case of odd number of the eigenvalues 

So far, we considered the case of an even number + 1 of eigenvalues Xs, s = 0,1, . . . , N 

with N — 1, 3, 5 Using the Christoffel transform, it is possible to obtain the 

para-Krawtchouk polynomials corresponding to an odd number of eigenvalues Xg, s — 
0,1, . . . , N — 1. Indeed, consider the polynomials 

X — Xn 

where A^ are given by (3.13). It is not difficult to check that 

^ ^ -Pn+i(a^jv) 

Pn{xN) 

and hence formula (6.1) defines new orthogonal polynomials Pn{x) obtained by the 
Christoffel transform [13] from the polynomials P„(x). The polynomials Pn{x) satisfy 
the recurrence relation 

Pn+l{x) + BnPn{x) + UnPn-l{x) = xPn{x) (6.2) 

with the coefficients [17] 

A 

Un = Un-l , Bn = -Bn+1 + " ^n- (6.3) 

A simple calculation yields 

~ _ n(N-n)((2n-Nr-{^-ir) 

{2n - Nf ^ ^ ^ 



and 



+ 7 7 - 1 iV , . . , , , 

Bn = - 1 + — — • 6.5 



The Jacobi matrix J with the recurrence coefficients Un, -Bn has size N x N and its 
eigenvalues Xg coincide with the eigenvalues of the matrix J apart from the final value 
xpf, i.e. 

Xs^s + ^{-f-l){l-i-iy). s^O,l,...,N-l, (6.6) 

Thus the polynomials Pn{x) defined by (6.2), (6.4), (6.5) have properties similar to those 
of the para-Krawtchouk polynomials. Clearly, the Jacobi matrix J generates a spin chain 
with PST property under the restriction (3.19). Indeed, Christoffel transforms hke (6.1) 
preserve the PST property and can be used to construct new examples of spin chains 
with PST (see [17] for details). 



Para-Krawtchouk polynomials 

Acknowledgments 



10 



AZ thanks Centre de Recherches Mathematiques (Universite de Montreal) for 
hospitality. The authors would like to thank M. Christandl and V. Spiridonov for 
stimulating discussions. The research of LV is supported in part by a research grant 
from the Natural Sciences and Engineering Research Council (NSERC) of Canada. 



Para-Krawtchouk polynomials 



11 



[1] C.Albanese, M.Christandl, N.Datta, A.Ekert, Mirror inversion of quantum states in linear 

registers, Phys. Rev. Lett. 93 (2004), 230502. 
[2] T. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, NY, 1978. 
[3] Ya.I.Granovskii, I.M.Lutzenko and A.S.Zhedanov, Mutual integrability, quadratic algebras and 

dynamical symmetry, Ann.Phys. (USA) 217 (1992), 1-20. 
[4] A.Kay, A Review of Perfect State Transfer and its Application as a Constructive Tool, Int. J. 

Quantum Inf. 8 (2010), 641-676; arXiv:0903.4274. 
[5] R. Koekoek,P. Lesky, R. Swarttouw, Hypergeometric Orthogonal Polynomials and Their Q- 

analogues. Springer- Verlag, 2010. 
[6] P. Letourneau and L. Vinet, Superintegrable Systems: Polynomial Algebras and Quasi-Exactly 

Solvable Hamiltonians, Ann. Phys. 243 (1995), 144-168. 
[7] A.F. Nikiforov, S.K. Suslov, and V.B. Uvarov, Classical Orthogonal Polynomials of a Discrete 

Variable, Springer, Berhn, 1991. 
[8] T. Shi, Y.Li , A. Song, C.P.Sun, Quantum-state transfer via the ferromagnetic chain in a spatially 

modulated field, Phys. Rev. A 71 (2005), 032309, 5 pages, quant-ph/0408152. 
[9] C.Smet and W.Van Asschc, Orthogonal polynomials on a bi-la,ttice, arXiv:1101.1817 math. 
[10] N.Stoilova, J. Van der Jeugt, An exactly solvable spin chain related to Hahn polynomials, SIGMA 

7 (2011), 033, 13 pages. 

[11] M. Rosenblum. Generalized Hermite Polynomials and the Bose-like Oscillator Calculus, in: Oper. 

Theory Adv. AppL, vol. 73, Birkhauser, Basel, 1994, pp. 369 396. ArXiv:math/9307224. 
[12] V.Spiridonov, L. Vinet, A.Zhedanov, Periodic reduction of the factorization chain and the Hahn 

polynomials, J. Phys. A: Math. Gen. 27 (1994), L669-L675. 
[13] G. Szcgo, Orthogonal Polynomials, fourth edition, AMS, 1975. 

[14] S.Tsujimoto, L. Vinet and A.Zhedanov, Dunkl shift operators and Bannai-Ito polynomials, 
arXiv:1106.3512. 

[15] S.Tsujimoto, L. Vinet and A.Zhedanov, From slq{2) to a Parabosonic Hopf Algebra, SIGMA 7 

(2011), 093, 13 pages; arXiv:1108.1603. 
[16] S.Tsujimoto, L. Vinet and A.Zhedanov, Dual -1 Hahn polynomials: "classical" polynomials beyond 

the Leonard duality, arXiv:1108.0132 (to be published in Proc.Amer.Math.Soc.) 

[17] L. Vinet and A.Zhedanov, How to construct spin chains with perfect state transfer, arXiv:1110.6474. 
[18] A.Zhedanov, The "Higgs algebra" as a quantum deformation of su{2), Mod. Phys. Let. A7 (1992), 
507-512. 



